Manipulating Bose-Einstein condensed atoms in toroidal traps 
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We consider Bose-Einstein condensed atoms confined in a toroidal trap. We demonstrate that 
under conditions of one-dimensional behavior, the density distribution of the atoms may be expo- 
nentially localized/delocalized, even for very small variations in the trapping potential along the 
torus. This observation allows one to control the atom density externally via slight modifications 
of the trapping potential. For similar reasons, small irregularities of the trap may also have a very 
pronounced effect on the density of the cloud. 
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I. INTRODUCTION 

One of the many interesting consequences of recent 
experimental advances in the field of cold atoms is the 
fact that it is now possible to realize states of reduced 
dimensionality. More precisely, so long as the energy 
of the atoms due to their interactions is much smaller 
than the spacing of the energy levels associated with mo- 
tion along some symmetry axis of the trap, the atomic 
state is dominated by the lowest state of the trapping po- 
tential, and this degree of freedom is effectively frozen. 
Thus, oblate/elongated traps can be used to achieve 
quasi- two/quasi-one dimensional behavior. As reported 
in Ref . Q , Ketterle and co-workers have achieved quasi- 
one-dimensional conditions. Using a variety of sophisti- 
cated techniques, experimentalists have produced traps 
spanning a wide range of other geometries, as reported in, 
e.g., Refs. 0, IE 0, |ETE • Several experimental and the- 
oretical groups have also studied quasi-one-dimensional 
Bose-Einstein condensates in magnetic waveguides |8| . In 
a detailed study, Leboeuf and Pavloff have considered 
the effects of perturbing obstacles on the propagation of 
Bose-Einstein condensed atoms through such magnetic 
waveguides @|. 

In the present study we focus on tight toroidal trapping 
potentials ,which have been realized in various laborato- 
ries 0, IE IE ■ Such traps can be used to create one- 
dimensional systems satisfying periodic boundary condi- 
tions. Since such systems have been under investigation 
for decades, it is clear that they are interesting from a 
theoretical point of view. Further, they may have numer- 
ous technological applications. 

In a one-dimensional toroidal trap, the atomic den- 
sity is homogeneous if the effective interaction between 
the atoms is repulsive |ldj |. It has been shown, however, 
that for a sufficiently strong attractive interaction, the 
atoms undergo a second-order quantum phase transition 
and form a localized density distribution 0, ^| . (This 
stands in contrast to a three-dimensional gas, which sim- 
ply collapses.) In the inhomogeneous phase the atoms 
benefit from the formation of a localized "blob" , which 
lowers their interaction energy. 

The issue to be considered here is the effect of a con- 



trolled change in the trapping potential V(9) along the 
torus as a consequence of small transverse irregularities. 
We will see that even small irregularities can materially 
alter this simple picture of localization. The effects to 
be considered were first encountered in early studies of 
radar. Hot spots were observed wherever waveguides 
were bent. It was soon recognized that any such non- 
uniformity in a quasi one-dimensional waveguide neces- 
sarily leads to exponentially localized, subthreshold res- 
onances (i.e., bound states) with a correspondingly high 
local energy density. 

In Sec. II we describe the problem to be considered 
in greater detail, and we present our model in Sec. III. 
Section IV contains our results for the ground state and 
the excitation spectrum of the gas. Section V contains a 
discussion of our results and presents the conclusions of 
this study. 



II. EFFECT OF A DISTORTED TOROIDAL 
POTENTIAL 

One might imagine that small irregularities in the one- 
dimensional potential, V(9), do not have any substantial 
effect. The situation is very different in the quasi one- 
dimensional systems considered here. Consider first a 
free particle moving in an infinitely long waveguide of 
small transverse size. The wave function for this particle 
can be approximated as the product of the lowest trans- 
verse eigenfunction and a function of the distance along 
the waveguide. The associated eigenvalue of the trans- 
verse Hamiltonian provides an effective potential for the 
resulting one-dimensional motion. A localized broaden- 
ing (or, equivalently, a bend) of an otherwise uniform 
waveguide will thus result in an attractive effective poten- 
tial. It is familiar from elementary quantum mechanics 
that every purely attractive potential in one dimension 
has at least one exponentially bound state. The ubiquity 
of bound states in infinite waveguides has been investi- 
gated thoroughly in an elegant paper by Goldstone and 
Jaffe [3. 

The situation is extremely similar for the toroidal ge- 
ometry considered here. Construct V(0) using the same 
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ansatz for the wave function, and set the zero of the en- 
ergy as the maximum value of V(6). Assuming a con- 
stant longitudinal wave function, an elementary varia- 
tional calculation immediately shows that the ground 
state energy of this system is necessarily negative if V{9) 
is non-uniform. There will always be classically allowed 
and forbidden regions with the usual attendant exponen- 
tial enhancement or suppression of the ground state wave 
function. 

For zero or weak interactions, the order parameter 
for the atomic system is simply the eigenfunction of 
the lowest energy state, and it is exponentially local- 
ized/delocalized, too. Thus, even small irregularities in 
V{9) in such one-dimensional systems render instability 
to an inhomogeneous state inevitable, provided only that 
the interaction is not too strong. Contrary to the quan- 
tum phase transition induced by an attractive effective 
interaction in a uniform toroidal trap, irregularities re- 
quire a sufficiently strong and repulsive interaction for 
the density to be (nearly) homogeneous. Furthermore, 
while the density variation induced by an attractive in- 
teraction in the uniform case is sinusoidal close to the 
transition point, irregularities in V{9) result in an ex- 
ponential localization/delocalization, which is generally 
more pronounced. 



III. MODEL 

To demonstrate these effects quantitatively, we con- 
sider the following form of V(9): 
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V(9) = 
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for |0| > tt/10 
for |0| < tt/10, 



(1) 



which corresponds to a broadening of the torus over 10% 
of its circumference. We assume the usual contact po- 
tential for the atom-atom interaction, V5 n t(r — r') = 
UoS(r — r') with Uq = 4:-irh 2 a/M. Here, a is the scat- 
tering length for elastic atom-atom collisions and M is 
the atomic mass. The Hamiltonian of the quasi-one- 
dimcnsional system thus becomes 



H/N = 



h 2 r f) 2 
2MR 2 / w <90 2 v ' 



v& t (0)V(0)tf(0)d0 + 

* t (0)* t (0)*(0)*(0)d0. (2) 



2nh Na 
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Here N 3> 1 is the number of atoms. The prefactor of 
the last term, 2nh 2 Na/(MRS), is equal to nnUo, where 
n is the (three-dimensional) density, n — N/ (2irRS), R 
is the radius of the torus, and S its cross section, with 
< R. The three energy scales that appear in our 
problem, are then (i) the energy for motion along the 
torus h 2 1 (2M R 2 ) , (ii) the depth of the potential Vq, and 
(hi) the interaction energy between the atoms, nUo- 



FIG. 1: The density |*(6»)| 2 of Eq. @ for three values of the 
coupling constant 7 = 0.1, 1.0, and 5.0, and for V(8) given by 
Eq. Jy . The two vertical lines show the range of the potential 
V(6). 



For Vq = 0, one finds an instability towards the for- 
mation of a localized state induced by an attractive in- 
teraction when n\Uo\ ~ h 2 /(2MR 2 ). More precisely, 
the critical value of the ratio, 7, between the poten- 
tial energy nUo and the kinetic energy H 2 /(2MR 2 ), is 
7 = ANaR/S = -1/2 PIE2. Finally, the crossover be- 
tween the phases of localized and homogeneous density 
that we consider here is given roughly by the condition 
nU Q ~ V . 



IV. RESULTS 



A. Lowest state 

To determine the order parameter ^(0), we first solve 
the eigenvalue problem -ffoV'm = E m ip m , where Hq is 
the Hamiltonian of the non-interacting problem. The 
eigenstates ip m are parity eigenstates. The actual value 
of Vo that we choose is 50 /ir 2 in units of h 2 /(2MR 2 ), 
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and it has only one bound state with an energy Eq w 
3.503h 2 /(2MR 2 ). Except for this lowest (to = 0) eigen- 
value, the eigenvalues of all the excited states are quite 
close to those of the undistorted torus, i.e., 



V 



2MR? 



(3) 



Clearly, this expression becomes more accurate for larger 
values of to. 

Having solved the eigenvalue problem, we then expand 
the order parameter in the eigenstates ij) m , 



(4) 



and calculate the coefficients c m variationally. For each 
value of the coupling constant 7 considered, we minimize 
the expectation value of the Hamiltonian subject to the 
constraint ^\ \c m | 2 = 1 corresponding to a fixed number 
of atoms. 

Figure 1 shows ^(S)! 2 for three values of 7 = 0.1, 1.0, 
and 5.0, where we kept the lowest five states (of even 
parity - only even parity states contribute to Vf). The 
highest eigenstate, with quantum number too, that must 
be included in the expansion of the order parameter has 
to satisfy the condition E mo 3> uUq. This condition is 
comfortably met even for the largest interaction strength 
7 = 5.0 considered, since E 5 « 29.620?i 2 / (2 MR 2 ). 

For small values of the interaction energy, uUq <C Eq, 
the dominant component of the order parameter is the 
lowest eigenstate of the single-particle problem. This is 
an exponentially decaying state in the classically forbid- 
den regions of the torus, \9\ > ir/10. Clearly, a simi- 
lar result will be obtained for any non-uniform choice of 
V(6). As the coupling increases, the density distribu- 
tion becomes wider. For sufficiently strong interactions, 
nUo ^S> Vo, the density ^(S)! 2 becomes homogeneous, 
reaching the limiting value l/(27r). 



B. Excitation spectrum 

The excitation spectrum is also affected as one goes 
from the limit of zero interactions (with the order pa- 
rameter being a localized state) to very strong interac- 
tions (with the order parameter being a homogeneous 
state). For weak interactions, the excitations are single- 
particle like. Atoms are excited from the eigenstate of the 
non-interacting problem with quantum number to = to 
some excited state with quantum number to =/= 0. In the 
absence of interactions the corresponding excitation en- 
ergies are approximately equal to the quadratic law of 
Eq. 0. 

For weak interactions (7^1) the excitation spectrum 
can be calculated using perturbation theory. The total 
energy of the lowest state is 
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FIG. 2: The points (that are connected with straight, dashed 
lines) show the excitation spectrum, Eq. Q, for 7 = 0.1 for 
the lowest five excited states of even and odd parity. The 
solid lines shows the free-particle excitation energy E m — Eq. 
The energy is measured in units of Ti 2 /(2MR 2 ). 



where Iijki = J ip*4'j' l Pk?Pi d0. Exciting one atom to the 
state with quantum number to, the total energy becomes 
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(N - 1)Eq + E m + 

[(N - 1)(N - 2) 7 0000 + 4(N - 1) I QmQri 
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The factor of four in the last term comes from the di- 
rect and the exchange terms of the interaction energy 
for identical bosons. From Eqs. JSJ and JfJJl we get for 
AE m = E m - Eq that 



AE m = E m - E + 27T7 
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Figure 2 shows AE m for the five lowest excited eigenstates 
of the non-interacting system of even and odd parity, for 
7 = 0.1. The excitation energies of the first two excited 
states (i.e., the even and odd parity states with \m\ = 1) 
are » 1.907^ 2 /(2Mi? 2 ) and w 2.443^ 2 /(2Mi? 2 ), respec- 
tively. Clearly for large enough \m\ the excitation spec- 
trum is quadratic. On the other hand, for the low-lying 
excited states there are deviations from the quadratic de- 
pendence due to the interactions, and also due to the de- 
viations of the eigenenergies from the formula of Eq. , 
which are more pronounced for low to. 
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The opposite limit of strong interactions, 7 3> 1, is 
more straightforward. Even in the presence of some V{&) 
the density tends towards homogeneity in this limit (as 
shown in the top panel of Fig. 1). The excitation spec- 
trum of the homogeneous state, | V I / (^)| 2 = l/(2n), is 

PUS 

For 7 ^> 1 the above formula implies that A£ m scales as 
I m I in our problem, 

A^ = ^v^H- 0) 

The above equation implies a speed of sound c = 
(R/h) dA£ m /dm = hy/Wy/(2MR), or Mc 2 = nU . This 
result can also be obtained from the energy per unit 
length, e(er), according to the formula Mc 2 = ad 2 e/dtj 2 , 
where a = N/ (2irR) is the atom density per unit length. 
In the limit of strong interactions, where \^f\ 2 = l/(2w), 
the quadratic contribution to e(cr) is given by the in- 
teraction energy, i.e., the last term in the Hamilto- 
nian H of Eq. ©, e(cr) = 2irh 2 ao 2 / (M S) . Thus, c = 
hy/2~7/(2MR), in agreement with Eq. ©. 

V. DISCUSSION AND CONCLUSIONS 

Depending on the various parameters and on the ex- 
perimental demands in constructing axially-symmetric 
toroidal traps, it is possible to explore relevant parts of 
the phase diagram of quasi one-dimensional atomic gases. 
We have pointed out, however, that the ground state of 
an atom in a toroidal trap will always have some degree 



of exponential localization whenever the trapping poten- 
tial deviates from perfect uniformity. Although we have 
illustrated this fact with one specific form of V(9), the 
result is general. Even apparently small modulations in 
V(9) can give rise to the localization demonstrated here 
and have the potential to obscure the quantum phase 
transition described above for an effective attractive in- 
teraction between the atoms. Furthermore, the excita- 
tion spectrum is also affected by the presence of some 
variation in V{6). More specifically, this is free-particle- 
like for zero or sufficiently weak interactions, and it is 
phonon-like for sufficiently strong interactions. 

Our study leads to three conclusions: (i) In realistic 
toroidal traps, even weak irregularities can result in a 
density distribution that is exponentially localized. Care 
must therefore be exercised in comparing experimental 
results with theoretical predictions assuming a uniform 
trapping potential, (ii) The excitation spectrum (and 
thus the dynamics) of the gas is also affected by varia- 
tions in V(9), and caution is again called for. (iii) The 
controlled modification of the toroidal potential and/or 
tuning of the coupling constant of the atomic interaction 
permit engineering of the shape of the atom density in 
a rather dramatic way. For example, the creation in the 
torus of a potential resembling our V(9) (e.g., using a 
laser) can change the density distribution from almost 
homogeneous to highly localized. Since such modifica- 
tions of toroidal traps can be made with relative ease 
and that the coupling between the atoms is easily tun- 
able via Feshbach resonances, our results may have useful 
applications. 

GMK acknowledges financial support from the Euro- 
pean Community project ULTRA- ID (NMP4-CT-2003- 
505457). 
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